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This is survey of some recent results connecting random matrices, non- 
colliding processes and queues. 



1 Introduction 

•5 

It was recently discovered by Baik, Deift and Johansson || that the asymptotic dis- 
tribution of the length of the longest increasing subsequence in a permutation chosen 
uniformly at random from S n , properly centred and normalised, is the same as the 
asymptotic distribution of the largest eigenvalue of an n x n GUE random matrix, 
| properly centred and normalised, as n — > oo. This distribution had earlier been iden- 



tified by Tracy and Widom IB3 in the random matrix context, and it is now known 
q ■ as the Tracy- Widom law. 

CN ■ 

O ■ The length of the longest increasing subsequence in a random permutation can be 

thought of as a 'last-passage time' for a certain directed percolation problem on the 
plane; this directed percolation problem is closely related to the one-dimensional 
totally asymmetric exclusion process (with low density of particles travelling at high 



speed) or equivalently, an infinite series of M/M/l queues in tandem (with low density 
of customers and high service rates). On the other hand, the eigenvalues of a GUE 
random matrix of dimension n have the same law as the positions, after a unit length 
of time, of n independent standard Brownian motions started from the origin and 
conditioned (in the sense of Doob) never to collide. These remarks are for the purpose 
of convincing the reader that there might be some interesting connections between 
random matrices, non-colliding processes and queues. Indeed there are, and that is 
the topic of this paper. 

Let us concentrate on the following, more exact, connection between directed percola- 
tion and random matrices which was more recently observed by Baryshnikov || and 
Gravner, Tracy and Widom ||26|| . Let B = (B\, . . . , B n ) be a standard n-dimensional 
Brownian motion, and write B(s,t) = B(t) — B(s) for s < t. 

Theorem 1 [Baryshnikov; Gravner- Tracy- Widom] The random variable 

M n = sup {B n (0,si) + ..- + Bi(s n _i,l)}, (1) 

0<«i<—<«„_i<l 

has the same law as the largest eigenvalue of an n-dimensional GUE random matrix. 



The proofs of this fact given in || and [26| are based on the Robinson-Shensted-Knuth 
(RSK) correspondence, and do not make use of the facts that 



(a) M n has a queueing interpretation [|4j, and 



(b) the largest eigenvalue of a GUE random matrix has an interpretation in terms 



of non-colliding Brownian motions [[H], |25 



In ||47||, a proof is given of a more general result which is based entirely on these 



interpretations. This more general result identifies a path-transformation T n (B) of 
B, which has the same law as that of n independent Brownian motions conditioned 
(in the sense of Doob) never to collide. This process, which we denote by B, is the 
eigenvalue process associated with Hermitian Brownian motion and B(l) has the same 
distribution as the eigenvalues of an n-dimensional GUE random matrix pi, 25 1 - The 



transformed process T n (B) has the property that its largest component at time 1 is 
given by M n , and so Theorem [l] follows. 

At the heart of the proof, which will be outlined in this paper, is a celebrated theorem 
of classical queueing theory, which states that, in equilibrium, the output of a stable 
M/M/l queue is Poisson. This result is usually attributed to Burke, who gave the 
first proof in 1956. In what follows, we shall also refer to it as an 'output theorem'. It 
follows from the reversibility of the M/M/l queue. There is an easy extension of this 
theorem, which can be proved by similar reversibility arguments; when phrased in the 
setting of 'max-plus algebra', this extension immediately yields the two-dimensional 
result. The result in higher dimensions is then obtained by considering a series of 
queues in tandem and applying an induction argument. The statement of Theorem [l] 
seems considerably less mysterious (to me at least!) in the setting of queues and 
non-colliding processes. 

In the case n = 2, the fact that T n (B) and B have the same law is equivalent to 
Pitman's representation for the three-dimensional Bessel process. In the case n = 3, 
it yields a representation for planar Brownian motion conditioned to stay forever in a 
wedge of angle 7r/3; a partial converse of this result was discovered earlier by Biane 0. 

As we remarked above, the process B has the same law as the eigenvalue process 
associated with Hermitian Brownian motion. Bougerol and Jeulin M have (inde- 
pendently) obtained a similar representation for this process, by completely different 
methods, which is also consistent with Theorem [TJ. Their results are presented in 
the more general setting of Brownian motion on symmetric spaces. The relationship 
between these two representations will be discussed elsewhere. 

There are certain symmetries in the max-plus algebra which play a crucial role, and 
these symmetries also exist in other algebras, including the conventional algebra on 
the reals. As a consequence, there are analogues of these output theorems in many 
different settings, and it seems that there are some general phenomena at work. In 
all cases, symmetry and reversibility play a key role. We will present some of these 
examples and their implications. 

We will also describe briefly how output theorems and the study of tandem systems 
can be used to obtain first order asymptotic results for various directed percolation 



and directed polymer models. 



The outline of the paper is as follows. In the next section we present some background 
material on random matrices and non-colliding Brownian motions. In section [3] we 



give a precise statement of the representation theorem obtained in [fi7 |, and briefly 



discuss the special cases n = 2 and n = 3. In section [|we describe the output theorem 
for the M/M/l queue, how it is proved using reversibility, and how it can be extended 
using similar reversibility arguments. In section [5] we show how the extended version 
of Burke's theorem can be used to obtain a representation for non-colliding Poisson 
processes, and how the representation for non-colliding Brownian motions follows. In 
section [| we discuss output theorems generally and give some examples. In section [?| 
we describe how output theorems can be used to obtain first order asymptotic results 
for various directed percolation and directed polymer models. 



2 Random matrices and non-colliding Brownian 
motions 

Recall that a n x n GUE random matrix A G <C nxn is constructed as follows: it 
is Hermitian, that is, A = A*{= (A) l )\ the entries {Aij, i < j} are independent; 
on the diagonal An are standard real normal random variables; below the diagonal, 
{Aij, i < j} are standard complex normal random variables, that is, the real and 
imaginary parts of A^ are independent centered real normal random variables, each 
with variance 1/2; above the diagonal we set = A\j. Here, z = x — iy denotes the 
complex conjugate of z = x + iy. The joint density (with respect to Lebesgue measure 
on M. n ) of the eigenvalues of A is given by Weyl's formula (see, for example, fl39|): 



WA 1; ... ,A n ) = Z- i n(A,-A J )nex P (-A l 2 /2), (2) 



where Z = (27r) n / 2 U"Zi J 1 

This distribution has an interpretation in terms of non-colliding Brownian motions. 
The (Vandermonde) function h defined by 

Kx) = Yl(xj ~ x i) (3) 

i<j 

is harmonic on R n , and moreover, is a strictly positive harmonic function for Brownian 
motion killed when it exits the Weyl chamber 

W = {x G R n : x x < x 2 < ■ ■ ■ < x n }. (4) 

For x G M n , let F x denote the law of B started at x; for x G W, let P x denote the law 
of the /i-transform of B started at x, where h is given by (^). The laws F x and F x are 
related as follows. If T denotes the first exit time of B from W, and T% the natural 
filtration of B, then for A G J-"t, 

K(A) =E X (^j^l An{T>t} y (5) 



It is well-known, and easy to check using the Karlin-MacGregor formula (see, for 
example, |]35| ) that 



lim F x (X t G dy) = C t h(y) 2 ¥ (X t G dy), 



(6) 



where 



n-l 



n(n-l)/2 



Tlfi 

3=1 



(7) 



is a normalisation constant. Note that the RHS of (|6|) is equal to Jgue (y) when t = 1. 

There is a more sophisticated connection at the process level. Note that we can define 
the law P 0+ , for A G T t = a(B u , u>t),t>0, by 



(8) 



where 9 is the shift operator (so that t A G T ). The fact that this is well-defined 
follows from (^). Now, Hermitian Brownian motion is constructed in the same way 
as a GUE random matrix, but with Brownian motions instead of normal random 
variables. It is a fact that the law of process of eigenvalues of Hermitian Brownian 
motion is given by F 0+ (see, for example, pTl ER |3 



For related work on non-colliding diffusions and random matrices, see || [L4], 0, [35|] 
and references therein. 



3 A generalisation of Theorem J] 

Let £>o(M+) denote the space of cadlag paths / : R + — * R with /(0) = 0. For 
f,ge A)(R+), define f®gE D (R + ) and / g G A)(R+) by 

Cf® <?)(*) = in! l [f(s)+ g(t)-g(s)], (9) 

0<s<t 

and 

(fQg)(t)= svLv[f(s)+g(t)-g(s)}. (10) 

0<s<t 

Unless otherwise deleniated by parentheses, the default order of operations is from 
left to right; for example, when we write / ® g h, we mean (/ ® g) ® /i. Define a 
sequence of mappings : D (R+) fc — > A)(R+) /c by 

r 2 (/,</) = (/®<7, </©/), (ii) 

and, for > 2, 

r*(/i,... ,/*) = (/i®/ 2 ®---®/fc, 

r fc _i(/ 2 A, / 3 © (A ® /a), • • • , A © (A ® • • • ® A-i)))- (12) 



The operations <g> (and ©) arise naturally in a queueing context; they can also be 
regarded as operator products in the min-plus (resp. max-plus) algebra. We refer the 
reader to |2[ for more about max-plus algebra and its applications. 

Let B be a realisation of Po+, as defined in the previous section. Then the largest 
eigenvalue of a n x n GUE random matrix has the same law as B n (l), and Theorem [T] 
states that M n and B n (l) have the same law. 

Now observe that, if T n (B)i denotes the first component in the n- dimensional process 
T n (B), then 

r„(fl)i(t) = (Si® ■■■®B n ){t) (13) 
inf {b$ , + --- + B\ n) A: (14) 

0<si<-<s„_i<t I («n-l.*)/' V > 

thus, by symmetry, Theorem [l] states that T n (B)i(l) has the same law as the smallest 
eigenvalue of an n- dimensional GUE random matrix or, equivalently, the random 
variable -Bi(l). We remark that, although it is not immediately obvious, it can also 
be shown that 



T n (B) n = B n Q---QB u (15) 

from which it also follows that M n = T n (B) n (l) has the same law as B n (l). In 
fact E7|: 



Theorem 2 The processes T n (B) and B have the same law. 



The proof of this identity in law, given in |47j , is based on an extension of the output 



theorem for a stable M/M/l queue in equilibrium, and certain related symmetries in 
the max-plus algebra. 

In the next section we discuss the output theorem for M/M/l queues. In Section [| 
we show how one can deduce the analogue of Theorem |2] for non-colliding Poisson 
processes. This is first done in the case of unequal rates, where the non-collision 
conditioning is non-singular. The case of equal rates follows, and it is then purely a 
technical exercise to apply Donsker's theorem and recover the Brownian version. We 
remark that the Poisson case with equal rates is interesting in its own right, as it is 
closely related to the Charlier orthogonal polynomial ensemble. 

As we remarked in the introduction, the case n = 2 is equivalent to Pitman's represen- 
tation for the three-dimensional Bessel process, which states that, if X is a standard 
one-dimensional Brownian motion and M(t) = max 0<s<( X(t), then 2M — X is a 
three-dimensional Bessel process. The three-dimensional Bessel process, being the 
/i-transform of Brownian motion on K + with h(x) = x, can also be interpreted as 
Brownian motion conditioned to stay positive. To see the connection, first note that 
in this case, R = {B 2 — Bi)/ \[2 is a three-dimensional Bessel process and (B 2 +Bi) / v2 
is a standard Brownian motion, independent of R. On the other hand, 



^(f,9)2 + r 2 (f,g) 1 = g + f, 



and 

r 2 (/, 9)2 - r 2 (/, g)i = g / - / ® g = 2m - x, 
where x = g — f and m(t) = max 0<s<t x(s). 

We now consider the statement of Theorem |2| in the case n — 3. In this case Brownian 
motion in the Weyl chamber 

W 3 = {x e M 3 : Xi < x 2 < x 3 } 

can be mapped onto Brownian motion in a cone of angle 7r/3 on the plane. Let 

C = {z e C : |argz| < vr/6}, 

and define (ft : W3 — >• C by 

(x z -x x 2x 2 -x 1 -x 3 \ 

If B is a standard Brownian motion in M 3 , then is a standard Brownian motion 
in C, and moroever, (p(W 3 ) = C (see, for example, || fiS]). It therefore follows from 
Theorem || that 4>(T 3 (B)) is a Brownian motion conditioned (in the sense of Doob) 
to stay in C forever. As remarked by Biane H] , there is only one way to condition on 
this event. 

In fact, Biane proved that, if X is a standard Brownian motion in C conditioned 
to stay in C forever, and J t = inf u>t X(u), then 3(3 J — X) is a standard one- 
dimensional Brownian motion. Thus, if we set X = <f)(T 3 (B)), then S(3J — X) 
is a standard one-dimensional Brownian motion. The expression 9(3 J — X) can be 
simplified considerably, but I do not see how to write it as a linear function of B; thus, 
it is not clear at this point how to recover Biane's result directly from Theorem [2| 



4 An extension of Burke's theorem 

The stationary M/M/l queue can be constructed as follows. Let A and S be inde- 
pendent Poisson processes on M. with respective intensities < A < /i. For intervals 
/, open, half-open or closed, we will denote by A(I) the measure of / with respect 
to gL4; for I = (0,t] we will simply write A(t), with the convention that A(0) = 0. 
Similarly for S and any other point process we introduce. For t £ R, set 

Q(t) = sup{A(s,t]- S(s,t]] + , (16) 

s<t 

and for s < t, 

D(s,t}= A(s,t] + Q(s)-Q(t). (17) 

In the language of queueing theory, A is the arrivals process, S is the service process, 
Q is the queue-length process, and D is the departure process. With this construction 



it is also natural (and indeed very important for what follows) to define the unused 
service process by 



U = S-D. (18) 

We will use the following notation for reversed processes. For a point process X, 
the reversed process X is defined by X(s,t) = X(—t, —s). The reversed queue-length 
process Q is defined to be the right-continuous modification of {Q(— t), t G M}. 

Burke's theorem states that D is a homogeneous Poisson process with intensity A. On 
a historical note, this fact was anticipated by O'Brien |E| and Morse |EJ, and proved 
in 1956 by Burke 



m 



Independently, it was also proved by Cohen |]L5|| . In 1957, 
Reich |)(J gave the following very elegant proof which uses reversibility. The process 
Q is reversible (in fact, all stationary birth and death processes are reversible). It 
follows that the joint law of A and D is the same as the joint law of D and A. In 
particular, D, and hence D, is a Poisson process with intensity A. 

Burke also proved that, for each t, {D(s,t], s < t} is independent of Q(t). This 
property is now called quasi-reversibility. Note that it also follows from Reich's re- 
versibility argument. Discussions on Burke's theorem and related material can be 
found in the books of Bremaud || [K| , Kelly and Robert |)1 . 



Set 



T = A + U. 



(19) 



Theorem 3 The point processes D and T are independent Poisson processes with 
respective intensities A and [i. 

Proof. First note that, given Q, U is a homogeneous Poisson process with intensity 
fi on the set / = {s G K. : Q(s) = 0}, and if we let V be another Poisson process with 
intensity fi on the complement of /, which is conditionally independent of U given Q, 
then (unconditionally) N = U + V is a homogeneous Poisson process with intensity 
/ionl which is independent of Q. Now, (A, S) can be written as a simple function 
of (Q,N), (A, S) = (p(Q,N) say. By construction, we have (D,T) = (p(Q,N). Now 
we use the reversibility of Q and iV to deduce that (D, T), and hence (D, T), has the 
same law as (A, S), as required. □ 

Note that 

Q(t) = sup[D(t, u) - T(t, u)}. (20) 

u>t 

We also have, on {Q{0) = 0}, 

{(£>(*), T(*)), t>0} = {T 2 (A,S)(t), t>0}. (21) 

Theorem |3] has the following multi-dimensional extension, which relates to a sequence 
of M/M/l queues in tandem. Let A, Si, . . . , S n be independent Poisson processes 



with respective intensities A, //i, . . . ,/i n , and assume that A < minj< n /ij. Set D = A 
and, for k > 1, t G IR, set 



and for s < t, 



Q fc (t) = supfDfc-^s, t] - S*(s, *]]+ (22) 



£> fc (s, t] = £> fc _i(s, t] + g fc (s) - Q fc (t), (23) 
T fc (s,t] = S , fc (s,t]-Q fc (s) + g fc (t). (24) 



Theorem 4 T7ie pomt processes D n ,Ti,... ,T n are independent Poisson processes 
with respective intensities A, fix, . . . ,[i n . 

Proof. By Theorem ^, D\,T\ and S2 are independent Poisson processes with respec- 
tive intensities \, fi± and fi 2 . Applying Theorem |3| again we see that D 2 and T 2 are 
independent Poisson processes with respective intensities A and fi 2 , and since D 2 and 
T 2 are determined by D\ and S2 they are independent of T±. Thus D 2 ,Ti,T 2 and S3 
are independent Poisson processes with respective intensities A,/ii,/x 2 and /X3. And so 
on. The condition A < minj<„ //, ensures that this procedure is well-defined. □ 



By repeated iteration of (£2|) and fl23|) , we obtain (almost surely) 

Qi(0) + • ■ ■ + Q n (0) = sup[A{s) - (5i ® • • • ® S„)(s)]. (25) 

s>0 

Remark. This formula, and variants of it, has been known for some time now. The 
first formula of this kind was observed by [[IT] and later extended in ^3] . Formulas 



of this kind can be used to compute first order asymptotics for directed percolation 
(and polymer) variables, following a program introduced by Seppalainen |52| (see |P[ 



for a survey). See also \ Tf[ f4"6fl for applications in a Brownian context. We give a brief 



outline of these ideas in Section [7|. 
Iterating (p3|) we obtain, for each k < n, 

D k {t) + Q x {t) + ■■■ + Q k {t) = A{t) + Qx(0) + • • • + Qfc(O). (26) 
We also have, by (p0|), 

Q fc (i) = sup[L> fe (t, u) - T fe (t, u)]. (27) 

u>t 

Applying this repeatedly we obtain 

Qi(O) + • • • + Q n (0) = sup[Z> n (t) - (Ti ® • • • ® T„)(t)]. (28) 

Note that, on {Qi(0) + • ■ ■ + Q n (0) = 0}, 

D n (t) = (A® Si® •••®S r «)(t), (29) 

and 

r fc (t) = (S A (A® S x ® • • • ® S*_i))(t), (30) 

for i > 0, k < n. 



5 Representations for non-colliding processes 



We will start by showing how the two-dimensional analogue of Theorem [| for Poisson 
processes is an immediate consequence of Theorem the extension of Burke's theorem 
given in the previous section, and the symmetry formula fl2~0|). As before, let A and 
S be Poisson processes with respective intensities A and /i. 



Theorem 5 The conditional law of {(A(t) , S (t)) , t > 0} given that A(t) < S(t) for 
allt>0 is the same as the unconditional law of {T 2 (A, S)(t), t > 0}. 



Proof. By Theorem ^ and the formula (|20|), the conditional law of {(A(t), S(t)), t > 
0} given that A(t) < S(t) for all t > is the same as the conditional law of 
{(D(t),T(t)), t > 0} given that Q(0) = 0. But when Q(0) = 0, (D(t),T(t)) = 
T 2 (A, S)(t) for t > 0. Moreover, by (|l^) and the independence of increments of A 
and 5*, {T 2 (A, S){t), t > 0} is independent of Q(0). Therefore, the conditional law 
of {(A(t),S(t)), t > 0} given that A(t) < S(t) for all t > is the same as the 
unconditional law of {T2(A S){t), t > 0}, as required. □ 

Note that we can immediately deduce a discrete version of Pitman's theorem, as 
follows. If we set X = S — A, and M(t) = max < s <( X(s), then X is a simple random 
walk with positive drift, and the law of 

2M-X = T 2 (A, S)z-r 2 (A, S)i 

is the same as that of X conditioned to stay non-negative. 

Now let = (Ni , . . . , Ni^) be the counting functions of n independent Poisson 
processes on M + with respective intensities fii < fi 2 < • • • < That is, Njf k \t) 
is the measure induced by the k th Poisson process on the interval (0, t], with the 
convention that Njf k \o) = 0. The following result was obtained in [f47| . 

Theorem 6 The conditional law of given that 

N^ l] (t) <■■■< N^\t), for all t > 0, 
is the same as the unconditional law o/r n (A r ^). 



Proof. We prove this by induction on n; the proof for n = 2 is given above. As- 
sume that Theorem ^ is true as stated for a particular value of n, and moreover 
holds for any choice of /xi < ... < /i„. In the above setting we have, by Theo- 
rem f|, that D n , Ti, . . . ,T n are independent Poisson processes with respective intensi- 
ties X, fii, . . . , fin- Assume that X < fii < ■ ■ ■ < Hn- By the induction hypothesis, the 
conditional law of 

{(D n (t),Ti(t),... ,T n (t)), t>0}, (31) 
given that Ti(t) < ■ ■ ■ < T n (t) for all t > 0, is the same as the (unconditional) law of 



{(£>„(*), r n (T!,... ,T n )(t)), t>0}; 



(32) 



therefore, the conditional law of 



{(Dnit),^),... ,T n (t)), t>0}, (33) 

given that D n {t) < T\(t) < ■ • ■ < T n (t) for all t > 0, is the same as the conditional 
law of 

{(D n (t),T n (T u ... ,T n )(t)), t>0}, (34) 

given that D n {t) < (7\ ® • • • <g> T n ){t) for all t > 0. But, by this is precisely the 

condition that Qi(0) + ■ • • + Q n (0) = or, equivalently, Qi(0) = • • • = Q n (0) = 0, 
and in this case we have, by (^) and (|3~0D , 



(D n (t),T n (T u ... ,T n )(t)) = r n+1 (A,St,... ,S n )(t) (35) 

for t > 0; since this latter expression, by independence of increments, is independent 
of Qi(0) + h Q n (0), we are done. □ 

Now let N = (Ni, . . . , N n ) be a collection of independent unit-rate Poisson processes, 
with N(0) = (0, . . . ,0). The function h given by (|3|) is a strictly positive harmonic 
function for the restriction of the transition kernel of N to the discrete Weyl chamber 



E = W fl 7j n . This is well-known to harmonic analysists; a proof is given in ||36|| . 
(In fact, h is harmonic for any random walk, or any Levy process, with exchangeable 
increments, provided the required moments exist.) Let x* + N be a realisation of the 
corresponding Doob /i-transform, started at x* — (0, 1, ... , n— 1) G -E (so that iV(0) = 
(0, . . . ,0)). Apart from providing a convenient framework in which to apply Donsker's 
theorem and deduce Theorem from Theorem 0, the process N is interesting in its 



own right. In [36] it is shown that the random vector N(l) is distributed according 



to the Charlier ensemble, a discrete orthogonal polynomial ensemble. That is, 

P(N(t) = y) = C t h(x* + y) 2 P{N{t) = y), (36) 

where Ct is the same (!) normalisation constant as in the Brownian case, given by (0). 
Thus, the next result, which follows from Theorem 0, yields a representation for the 



Charlier ensemble. For more on discrete orthogonal polynomial ensembles, see ]32 . 
Theorem 7 The processes N and T n (N) have the same law. 

Some elementary potential theory is needed to deduce this from Theorem since 
the function h is not the only positive harmonic function for the restriction of the 
transition kernel of N to the discrete Weyl chamber E = WnZ/ 1 . We refer the reader 
to |47J for details of the proof, and to [^| for the necessary asymptotic analysis of 



the Green's function. We remark that, for any y e IR™, the function defined by 
g y (x) = Schur x (?/) is a positive harmonic function for iV on E, and note that g y is 
a multiple of h if all the components of y are equal. This follows from the analysis 
of the Green's function presented in [^] (alternatively it can be deduced from more 
general representation theoretic arguments, as in 0). 

Theorem ^| (the Brownian version) follows from Theorem [7] by careful application of 
Donsker's theorem (see [^7] for details). 



The analogues of Theorems ^| and |7] for discrete-time walks are presented in [36] ; in 
this case, replace 'Poisson' by 'binomial', 'Charlier' by 'Krawtchouk', and modify the 
definition of T n . The proofs are similar, but there are additional complications due 
to the fact that the walkers can jump simultaneously; the asymptotic analysis of the 
Green's function is also more difficult in this case. 



6 Output theorems generally 

Theorem seems to be a special case of a more general phenomenon. In this section 
we will present several examples of 'output theorems'; in all cases, reversibility and 
symmetry play a key role. 

The analogue of Theorem [^ holds for Brownian motions with drift. Suppose now that 
A and S are independent Brownian motions, indexed by R, with respective drifts 
A < fi. That is, for each s, {A(s + 1) — A(s), t > 0} is a Brownian motion with drift 
A. Set A(s,t) = A(t) - A(s), S(s,t) = S(t) - S(s), 

Q(t) = sup[A(s,t) -S{s,t)}, 

s<t 

D(t)=A(t)+Q(0)-Q(t), 
T(t)=S(t)-Q(0)+Q(t), 
for i 6 I. Then the pair (D,T) has the same law as (A, S). Several proofs of this 



fact are given in [|6|; an equivalent result was earlier presented in p9fl (see below). 



Note that Theorem || can also be stated as follows, by considering the process X = 
S — A and its transform X = T — D. 

Theorem 8 Let X be a simple random walk (in continuous time) with positive drift, 
indexed by R ; and set 

Q(t)= sup X(s)-X(t), 

— oo<s<t 

X{t)=X{t)+2[Q{t)-Q{0)\. 
Then X has the same law as X . 

Note that, in this setting, the 'symmetry formula' (20) becomes 



Q(t) = swpX(t)-X(u). (37) 

u>t 

Again, the Brownian version of this result holds: this was the Brownian version of 



Burke's theorem presented in [29 



The analogue of Theorem |] can also be shown to hold for Brownian motions with 
drifts, by exactly the same argument given in the proof of Theorem |j. 

We will now describe some output theorems, which are completely analogous to those 
described above, but in the usual algebra (as opposed to the max-plus algebra). These 



were given in [46 1, and motivated by the recent extensions of Pitman's 2M—X theorem 
obtained by Matsumoto and Yor |37| in the context of exponential functionals of 
Brownian motion. In words, the analogues of Theorems ^ and |8| hold for Brownian 
motions if we replace 'sup' by 'log f exp' in the definition of Q. More precisely, suppose 
that A and S are Brownian motions with respective drifts A < fi and we set, for iGR, 

Q(t) = log ! expL4(s, t) - S(s, t)]ds 

J — oo 

D(t) = A(t) + Q(0)-Q(t), 

T(t)=S(t)-Q(0)+Q(t). 

Then the pair (D, T) has the same law as (A, S). As discussed in [[PJ, symmetry and 
reversibility play a big role here. In particular, the process Q is reversible, and the 
analogue of the symmetry formula ( |5T| ) holds: 



Q(t) = log / exp[D(t, u) - T(t, u)}du. 



(3? 



Let us return to the max-plus version for a moment: intuitively, the statement of 
Theorem |3| can almost be seen as a direct consequence of the reversibility of the 
process Q combined with the symmetry formula (|20|); the only thing which makes this 
difficult is the non-invertibility of the transformation. This problem does not exist in 
the usual algebra, and in fact the statement that {D,T) has the same law as {A,S) 
(when everything is defined using log J exp) is actually equivalent to the statement 
that Q is reversible. For a more precise statement of this equivalence, which holds 

The version of Pitman's 2M — X theorem obtained in this 



quite generally, see 



context by Matsumoto and Yor j37| states that, if X is a Brownian motion with drift 
and 

2M(t) = log f e 2X(s) ds, 
Jo 

then 2M—X is a certain diffusion process. Note that by following the arguments given 
at the start of the last section to deduce Pitman's theorem for random walks from 
the output theorem [| we can (formally) recover the Matsumoto- Yor representation 
and, moreover, see that, in some suitable sense, the law of the diffusion 2M — X can 
be interpreted as the conditional law of X, given that 



1. 



Similarly, the analogue of Theorem |] can be shown to hold in the log J exp-world for 
Brownian motions with drifts, by exactly the same extension argument given in the 
previous section. Note also that the arguments given in the proof of Theorem |6] can 
also be (formally) carried over to obtain a log J exp-analogue of Theorem 0. This has 
the following form: define a mapping IT n analogously with T n but replacing sup by 
log J exp (and inf by — log J exp(— )), in the definitions of ® and 0. Then, by these 
formal arguments, the law of U n (B) can be interpreted as the conditional law of B, 
given that 

POD 

/ exp (s) - B i+1 (s)]ds = I, 
Jo 



for each i = 1, . . . , n — 1. 



In [p8 |, the following non-Markovian version of Pitman's theorem is obtained using 
similar ideas: reversibility arguments do not require the Markov property. 

Theorem 9 Let k > 0) be a Markov chain on {— 1,+1} with Co — 1 tran- 
sition probabilities P(^k+i = 1| £fc = 1) = a -P(6c+i = — 1| 6c — — 1) = b < a. Set 
Xq = 0, X n = £i + • • • + and M n = maxo<fc<„ X^. Tae process 2M — X nas iae 
same /aw as iaai o/X conditioned to stay non-negative. 

The continuous analogue of the above theorem, with £ replaced by a zig-zag process, 



is also given in [2S 



We end this section with one final example, quite different from the above, to give 
an indication of how generally these output theorems, and the corresponding rep- 
resentation theorems, hold. It is related to the first-order autoregressive process. 
Let (b n , n G Z) be a sequence of independent standard normal random variables, 
and fix < a < 1. Consider the stationary process (X n , n G Z) defined by 
= o a 3 ^n~j- Then X n+ \ — aX n + b n for all n. Now observe that, if we 

set 

b n = X n — aX n+ i, 

then X n = aX n+ i + b n and X n = Ylf=Q a ~^n+j- Now, by the symmetry of this 
construction, and the fact that the process X is reversible, we see that: 

Theorem 10 The sequences b and b have the same law; that is, (b n , n G Z) is a 
sequence of independent standard normal random variables. 



This is the output theorem. For the analogue of Pitman's theorem in this case, set 
Y (x) = x and, for n > 0, = aY^ x) + b n , 

Ax) _ y(x) _ y(x) 
u n ~ 1 n U1 n+1- 

Consider the conditional law of (b n , n > 0) given that £^L a Tl &n = (It is easy to 
see that this is well-defined.) From the output theorem, and the symmetry formula 
X = £^L a~ n b n , this is the same as the law of (b n , n > 0) given that X = x. But, 
when X = x, b n = dn and, moreover, (dn \ n > 0) is independent of X ; it follows 
that 

Theorem 11 The law of the sequence (d ( n\ n > 0) is the same as the conditional 
law of (b n , n > 0) given that £^1 a n b n = x. 



7 First-order asymptotics for directed percolation 
and directed polymers 

The formula 

Qi(0) + • • • + Q„(0) = sup[A(s) - (St ® • • • ® S n )(s)}, (39) 

s>0 



introduced in Section f|, is very useful for computing first-order asymptotic results for 
directed percolation variables. 



In Section [|, A, Si, . . . , S n are independent Poisson processes with respective rates 
A, /ii, . . . , fi n , with A < minj /Zj. In this case, the random variables Qi(0) are indepen- 
dent and geometrically distributed with respective parameters A//ij. Assume A < 1. 
If we let Hi = 1, for each i, and divide (^) by n, we can (in principle) let n — > oo to 
obtain: 

r- =sup[Aar-7(a;)], (40) 

1 — A x >0 

where, almost surely, 

j(x) = lim -(St ® • • • ® S n )(xn). (41) 

(Here we have implicitly used the reversibility of A, and of Si, for each i.) The exis- 
tence of 7 follows from Kingman's subadditive ergodic theorem. But (|40|) is essentially 
a Legendre transform. In particular, it can be inverted to obtain 

j(x) = (Vx- l) 2 l x >i. 



This ingenuous idea is due to Seppalainen |p2| , who used it to compute asymptotics 



for a certain last-passage Bernoulli directed percolation problem. See |44| for a survey 
on the application of this technique to a variety of directed percolation problems. In 
each context, there is an associated 'growth model', and these limit theorems can also 
be interpreted as determining the limiting shape of the associated growth model. See 
also [0TI for related work in this 



For the Brownian model, we have, almost surely, 

lim -(Si® ■■■®B n )(xn) = 2^/x, (42) 

rwoo n 

where Bi,B 2 ,... are independent standard Brownian motions. This is described 
in JiH] and proved with sharp uniform concentration estimates in [2~7|. The cor- 



responding corner growth model is also discussed in that paper. In this case, the 
formula (P^) holds with A,S\,... ,S n replaced by independent standard Brownian 
motions with respective drifts A, 1, . . . ,1 and the Qi(0) are i.i.d. exponentially dis- 
tributed with parameter A. Note that the limiting result (f|^) can also be deduced 
from Theorem [I] by applying standard results from random matrix theory. 

The formula (|39|) is also valid when 'sup' is replaced by 'log J exp'. In this context it 
is used in j|6| to compute (formally) the free energy density associated with a certain 



directed polymer in a random environment. More precisely, for (3 > 0, if 

Z n ((3) = / ds\ . . . ds n -i exp {(3(Bi(—n, s\) H h B n {s n - X , 0))} 

J — n<si<-<s n _i<0 

where Bi(s,t) = Bi(t) — Bi(s), then, almost surely, 

f{/3) = lim -logZ n (l3) = -g(-/3 2 ) - 2 log ft 

n^oo n 



where 

g(x) = sup[xy + #(?/)], 

y>0 

for rr < 0, and ^ = r'/r is the digamma function. 
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